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, ^ \ Phonon exchange with nuclei in the course of fusion reactions that occur in a sohd have not been 

analyzed previously. This problem has become of interest in connection with claims of observations 
of anomalies in metal deuterides. If the strong force interaction were dependent only on position 
D ' (and not spin or isospin), then the coupling with phonons can be developed directly. Since a 



nuclear interaction can change the lattice constituents, the initial and final state lattices can be 
different, and we must include this in the formulation. For more realistic strong force models with 
spin and isospin dependence, we can use correlated nuclear wavefunctions which are made up of 
. products of space, spin and isospin components. In this case, the spin and isospin algebra can 

■ be done analytically, producing channel-dependent potentials that are only space dependent. The 

Q ■ formulation that results can be used for quantitative estimates of phonon exchange. 
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I. INTRODUCTION 



Few claims have resulted in as much controversy as the claim of the observation of an 
excess heat effect of nuclear origin in metal deuterides made by Fleischmann, Pons, and 
Hawkins in 1989.- Reviews of work on this problem were conducted by the Department of 
Energy (in the United States) in 1989^ and in 2004.- The conclusions of the 1989 review 
were very negative. The 2004 review was quite different in how it was carried out, what 
material was reviewed, and in the results. On the specific question of the existence of an 
excess power effect, which was the central issue of the review material presented to DoE,- 
there were more positive comments than negative comments volunteered from the reviewers.- 
This is significant, as DoE did not charge the reviewers to consider this question specifically. 
Among many other comments that were made, several of the reviewers suggested that more 
work on this problem be submitted to mainstream journals. 

Our interest in this paper is not focused on the question of the existence of an excess 
heat effect, although this issue has motivated our investigations. We are interested instead 
in theoretical issues which must be better understood in order for theory to be more relevant 
to the problem. For most physicists, the theoretical problem was adequately addressed by 
Huizenga, who argued that three "miracles" were required for a theoretical explanation.- 

• The first "miracle" involves the question of reaction rate. The excess heat effect implies 
reaction rates on the order of 10^^ sec~^, tens of orders of magnitude larger than the 
fastest fusion rates possible in molecular D2 or HD. 

• The second "miracle" concerns the branching ratio. Fleischmann initially speculated 
that ^He was being made from two deuterons. Such a result is seemingly miraculous 
as the branching ratio for this reaction channel is lower by about seven orders of 
magnitude than the primary n+^He and p+t reaction channels. 

• The third "miracle" involves the requirement that the reaction energy be expressed 
through channels not involving energetic reaction products. 

If one believes that these issues are insurmountable, then one would require that any claim 
of excess heat of nuclear origin must be accompanied by measurable energetic reaction 
products. This view was adopted by Huizenga, and by the 1989 ERAB Panel. This view is 
also widely held among the physics community at this time. 

A. Motivation 

However, other views are possible. Fleischmann's initial speculation was that a new 
physical process was involved, one that behaves differently from what is described in the 
textbooks. It seems that if we are trying to understand whether an excess heat effect 
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can exist, we should be asking whether there can be interactions between nuclei and a 
condensed matter environment. For example, suppose that phonons are exchanged in the 
course of a nuclear reaction that occurs in a lattice. Could this phonon exchange result 
in a modification of the selection rules? And if phonon exchange occurs, would that not 
constitute a new channel into which reaction energy might go? Nuclear physics relies heavily 
on the notion of an equivalent vacuum reaction in order to understand quantitatively fusion 
reactions. Huizenga's three "miracles" refiect such a view. For example, the second "miracle" 
is essentially a statement of energy, linear momentum, and angular momentum conservation, 
all of which are conserved in vacuum. Is it obvious that conservation laws in a condensed 
matter environment always produce the same reaction pathways with the same branching 
ratios; and, if so, can this be proven? In the case of the third "miracle," the reaction energy 
must be expressed as an energetic product in vacuum, since no other kinds of channels exist. 
But is it obvious that this must also be the case in a condensed matter environment, where 
low energy channels do exist? Once again, can this be proven? In any event, all of these 
questions ought to be capable of being addressed by theory. 

The current lack of such a theory in the mainstream literature is underscored in comments 
made by some reviewers in the 2004 DoE review. For example, a reviewer with expertise 
in metal deuterides and phonons wrote: "To create a coupling between nuclear interaction 
and phonons at such a low energy region (namely, the electromagnetic interaction) is beyond 
one's imagination at the moment." Another reviewer with expertise in nuclear theory wrote: 
"I am convinced that simple order-of-magnitude estimates of this kind could quickly rule 
out any of the exotic mechanisms proposed..." This reviewer argues that the matrix element 
for the emission of a large number of phonons would scale as (kR)"^^, where k is the phonon 
wavenumber. Aside from the fact that the emission of a large number of phonons in a 
single process was not put forth as a proposed reaction mechanism in the review, the use 
(or misuse) by this reviewer of a phonon wavenumber in this way underscores the need for 
a relevant theory. We find ourselves then in a position where a better understanding of 
phonon exchange in fusion reactions is required for a sensible discussion of these and related 
issues. This motivates us to consider in this manuscript how to calculate phonon exchange in 
perhaps the simplest application; that of an interaction matrix element. Phonon exchange 
due to recoil in neutron capture reactions was analyzed by Lamb,- and this analysis has 
been used subsequently for phonon exchange in gamma reactions. In the Mossbauer effect, 
an anomalous zero-phonon exchange effect occurs if the recoil is sufficiently gentle.- In our 
matrix elements recoil is involved implicitly, and phonon exchange due to the local change 
in lattice structure will appear explicitly. 
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B. Excitation transfer scheme 

In the following sections of this manuscript, we develop a formulation suitable for includ- 
ing phonon exchange in nuclear reaction matrix elements. Once armed with this tool, we 
are in a position to develop calculations relevant to candidate reaction schemes that involve 
phonon exchange (of which quite a few have been proposed over the years). Our current 
interest is focused on a particular scheme in which coupling occurs between nuclei at differ- 
ent sites as a result of phonon exchange with a common phonon mode. In this scheme, two 
deuterons fuse locally with the reaction energy transferred elsewhere through a second-order 
off-resonant excitation transfer effect. The systematics of this effect in the idealized case 
of two-level systems coupled to a common oscillator is discussed in a recent analysis.- In 
essence, the reaction energy is transferred to excited states in host lattice nuclei through 
an initial slow excitation transfer step, and then subsequent rapid excitation transfer occurs 
among these nuclei. The reaction energy is transferred to the phonon mode in this scheme 
a few phonons at a time in a very large number of fast excitation transfer reactions. The 
scheme seems to be closely related to experiment, and the reaction rate predicted appears to 
be consistent with experiment to within uncertainties in the deuteron-deuteron screening. 
The evaluation of this scheme requires the ability to evaluate phonon exchange in fusion 
reactions as discussed in the present manuscript. 

Huizenga's arguments are based ultimately on fundamental notions relevant to vacuum 
reactions, local energy and momentum conservation, and the use of Golden Rule reaction 
rate estimates, as discussed above. In this new excitation transfer scheme, the energy of 
the coupled system is conserved, but local energy is transferred elsewhere. In addition, 
the associated reaction dynamics do not follow simple Golden Rule predictions. Huizenga's 
arguments simply do not apply to this kind of scheme. 

C. Overview of the paper 

This paper is organized as follows. In section II we discuss the basic lattice-nuclear 
coupling (assuming a simplified nuclear potential model that depends only on position), 
and argue that it is very similar to the more familiar electron-phonon and neutron-phonon 
couplings. However, realistic nuclear potentials depend on spin, isospin, and parity. In 
Section HI we show that for realistic nuclear potentials, by first doing the spin and isospin 
algebra analytically, the ideas developed in Section II are relevant. In Section IV, using the 
Hamada- Johnston nuclear potential, we explicitly carry out the program outlined in Section 
HI for one particular matrix element relevant to a deuteron-deuteron fusion reaction. This 
particular interaction is chosen because it is one of the simplest realistic nuclear potentials. 
In Section V we extend this vacuum result to include the lattice. We find in Section VI an 
explicit expression for one of the integrals that appears in the interaction matrix element in 
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terms of phonon coordinates. This result is closely related to the discussion in Section II. A 
summary and conclusions is given in Section VII. 

Although the essential issues involved in phonon exchange can be seen with the simplified 
position-dependent potential, we need to work with nuclear wavefunctions that depend on 
spin and isospin when using a more realistic nuclear potential model. The construction of 
these general four-body correlated nuclear wavefunction requires the use of certain aspects 
of representation theory, and may be unfamiliar to some of the readers. As a result, we have 
included material in a set of appendices that provide further discussion of the wavefunction 
construction and usage in matrix element calculations. 
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II. PHONON INTERACTIONS 



We are interested in the basic issue of phonon exchange in association with nuclear 
reactions in a lattice. The formulation that follows in subsequent sections is complicated, 
and it makes sense to review phonon exchange in more familiar circumstances. In condensed 
matter physics, phonon exchange is of interest in the study of electron scattering and neutron 
scattering. In both cases, phonon exchange comes about through a formulation in which 
the center of mass coordinates of the nuclei are taken to be phonon operators. 



A. Phonon exchange associated with electronic transitions 



For example, electron-phonon interactions can be developed ultimately starting from the 
Coulomb interaction written in the form^ 



(1) 



where the are the electron coordinates, and the Zj are the effective charges on the ions. 
The Rj are the nuclear center of mass coordinates, which are dynamical, and which we think 
of in terms of lattice operators {qm) 



Ri 



R-i(gi,g2, 



(2) 



In the event that we adopt a description of the lattice in terms of product states, then a 
matrix element of the Coulomb interaction might be written as 



M 



fi 



^?({9})^/({x}) 



5 IR-i Xfcl 



(3) 



In this way of thinking about the problem, we could define a phonon interaction by inte- 
grating over electronic coordinates 



I R7 — Xfe I 



^.({x}) 



(4) 



j,k 

Phonon exchange in this case can be developed then from lattice matrix elements of this 
operator 



Mf^ = (^?(W)l^)/.(W)l^f({g})> (5) 

In the event that the electronic transition results in a change in the local force constants, the 
initial and final state lattice may differ enough that this should be reflected in the matrix 
element. 
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B. Phonon exchange associated with neutron interactions 



A similar approach can be used in the case of neutron scattering. In the case of a simple 
position-dependent interaction 



(6) 



we can develop the phonon interaction according to 



VoJ2 5'(R.-rn) 




(7) 



In these equations r„ is the neutron coordinate. With a knowledge of the phonon interaction, 
we can analyze phonon exchange using Mfi matrix elements as calculated above. 



C. Phonon exchange in a simplified model for nuclear reactions 

We might choose to work with an approximate position-dependent nuclear interaction 
model that is a function of spatial coordinates alone 

J2M\^c.-rp\) (8) 

a</3 

where a and f3 refer to nucleon (neutron and proton) coordinates. The interaction matrix 
element in this case can be written as 



r/3|) 



(9) 



a<l3 

where the integration here is taken over individual nucleon coordinates. A nuclear reaction 
may result in the initial lattice being different than the final lattice. For example, if two 
deuterons reacted to form ^He-|-n, the initial lattice would contain the two deuterons while 
the final lattice would have ^He and a possibly free neutron. We note that in a normal 
version of this reaction, the final state products would fiy off as energetic nuclear products; 
we consider this possibility to be included within the formulation under discussion. 

However, this description does not yet make clear the situation with respect to phonon 
exchange. For this, we need to rewrite the integral in terms of phonon coordinates rather 
than in terms of nucleon coordinates. To make progress, we can express individual nucleon 
coordinates that are associated with nuclei in terms of the nuclear center of mass coordinates 
R and relative internal nuclear coordinates ^ 



*^({r}) = *^({R},{|}) 



(10) 



7 



Since a nuclear reaction may result in a change in nuclei between the initial states and the 
final states, the set of center of mass coordinates may be different. Consequently, the nuclear 
interaction matrix element can be written formally as 



M 



fi 



*^({Rn, {e})] * ^n(|r.-r^|)*f ({R^}, {C}) H {<-ri)d{e}d{e}dmd{Kf} 

(11) 



We integrate over all initial and final state center of mass and internal coordinates, and 
require that the individual nucleon coordinates be the same in the initial and final states. 

Once we have a description in terms of the nuclear center of mass coordinates, we can 
rewrite the integrations in terms of phonon coordinates. To simplify things, we might assume 
product states of the form 

*'^({R},u}) - ^"(W)*(m) (12) 

In this case, the matrix element becomes 



j [^'}{{qf})^s{W})] *Y.MK- ^p\)'^hWmm}) n ^'(^L - ^L)d{e}dW}d{q'}d{qf} 

(13) 

For a final lattice with ions, this can be thought of as 

Mf. = y"[vl>J(q/)]* ^;^,(q/,q^) vI/f(q^)5(3A^-3)(q/_A-q^-b)rfqW (14) 

in which the phonon interaction operator Vfi{q;^ , q*) is developed from an appropriate inte- 
gration over the internal nuclear coordinates. More simply, we may write 

Mf, = I [v[/f(q/)]* Vf,{qf,ct) ^f{q')dq' (15) 

together with the associated constraint 

q^ = Aq* + b (16) 

Due to the change in the lattice structure, the phonon mode structure can be altered.—*^ 
The associated matrix element written in terms of initial state and final state coordinates 
[equation f|T5|) ] is common in papers on Franck-Condon factors in polyatomic molecules.— 
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D. Discussion 

Phonon exchange in electron and neutron scattering can be implemented using the local 
nuclear center of mass operators as phonon operators. This leads to a straightforward 
description of matrix elements that include phonon exchange, as well as the development 
of phonon exchange operators. In the case of a nuclear reaction, we can define a nuclear 
interaction matrix element simply enough under the assumption of a position-dependent 
potential as long as we work with nucleon coordinates. However, when the integral is 
rewritten in terms of phonon mode amplitudes, the resulting expression is more complicated. 
The underlying ideas and approach are very similar. The important message here is that 
a nuclear interaction matrix element expressed in terms of the nucleons of the lattice is 
completely straightforward, and it is only the bookkeeping associated with the assembly of 
nucleons into nuclei, and nuclei into phonons, that adds some complication to the problem. 

There are basic differences between the new problem and the well known electron and 
neutron scattering problems which deserve to be noted. 

• The energy scales in the case of nuclear reactions are not matched; typical nuclear re- 
action energies are on MeV scale, whereas the maximum phonon energies are typically 
100 meV or less. 

• A nuclear reaction involving a change in the nucleon content of the product nuclei 
will likely produce phonon exchange, if for no other reason than the fact that such 
reactions will produce a change in the structure of the lattice. 
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III. FORMAL CALCULATION WITH A REALISTIC POTENTIAL 



For quantitative results, we will need to use a realistic nuclear potential. In atomic physics 
and solid state physics, we are familiar with wavefunctions that have both spin and space 
dependence. Nuclear wavefunctions have in addition isospin dependence. In the isospin 
scheme, a nucleon can be in an isospin down state | |) in which case it is a neutron, or it 
can be in an isospin up state | t) in which case it is a proton. It was recognized early on in 
nuclear physics that the nucleon-nucleon potential is of short range in space, and depends 
explicitly on both spin and isospin. Over the past 60 years, increasingly accurate nuclear 
potential models have been developed. Early models were based in part on field theoretical 
models for the one-pion exchange, and in part on few parameter empirical models that 
were fit to scattering data. In recent years, the most accurate nuclear potential models 
have been constructed from a diagrammatic analysis of an effective field theory.— For the 
purposes of the present paper, we will adopt an early potential model (the Hamada- Johnston 
potential^) , which has the advantage for our discussion of being relatively simple in form. 



A. Hamada-Johnston potential 



The Hamada-Johnston potential between nucleon 1 and nucleon 2 appears in the litera- 
ture written as 



Vhj{1,2) = Vc{l.2) + V^t(1, 2)^12 + V^l5(1, 2)(L.S) + Vll{1,2)L 



12 



(17) 



This potential is made up of four basic terms: a central potential Vc'-, a tensor interaction 
Vt"S'i2; a spin-orbit term Vls'(L.S); and a generalized centripetal potential term VllLi2- 
These are discussed further in Appendix A. The central and tensor terms have isospin 
dependence, whereas the spin-orbit and generalized centripetal terms do not. As an example, 
we consider the tensor term in more detail (which we use as an example in Appendix [D]) . 
Explicitly it is given by 



Vt{1,2)Si2 = {ti ■ T2) yrifLr 



12; 



^[o-i ■ ri2j (0-2 ■ ri2j 

O K (Ti ■ (T2 

^12 



(18) 



with fi = niT^c/h, where is the mass of the pion. The radial potential is parameterized 
according to 



0.08- 



, 3 3 \ e 
1 + - + - 



X 



e e 
1 + aT h &T 



-2x 



X 



x^ 



(19) 



where qt and Bt are parameters which have been fit separately for singlet-triplet spin and 
even-odd parity (angular momentum) channels. Hence, i/t has an additional implicit de- 
pendence on spin and parity: i/t = Vt where a can be even-triplet (et), even-singlet (es), 
odd-triplet (ot) or odd-singlet (os). 
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B. Wavefunctions 



Because realistic nuclear potentials have explicit dependences on spin and isospin, one 
requires nuclear wavefunctions that have well defined spin and isospin dependences. In the 
nuclear physics literature, there are a variety of approaches to this problem, ranging from 
determinantal wavefunctions built up from single nucleon orbitals to more complicated cor- 
related wavefunctions. For the development of the results in this paper, we will adopt a basic 
construction of correlated wavefunctions in terms of products of spatial, spin-dependent, and 
isospin-dependent pieces, each of which is determined by group theory. 

As is customary in nuclear physics, we choose our wavefunctions to satisfy the general- 
ized Pauli exclusion principle, which requires that the total wavefunction (including space, 
spin and isospin parts) has to be antisymmetric with respect to nucleon exchange. A group 
theoretical construction of the antisymmetric wavefunctions is possible by summing over 
products of many-particle spatial functions, spin functions, and isospin functions, each of 
which belong to appropriate representations of the symmetric group, S'(4).— In the case of 
the spin and isospin functions, the symmetric group representations correspond to eigen- 
functions of total spin and isospin, as a consequence of the Schur-Weyl duality.— Hence 
the antisymmetric wavefunctions have well-defined total spin and isospin. This is discussed 
in Appendix B. 

In general, such correlated wavefunctions can be written in the form 

^ = E^r^.(W)^.(M)^.(W) (20) 

i 

where Sj({cr}) are the spin functions, tj({r}) are the isospin functions, and 'ipjUr}) are the 
spatial functions; all of which belong to certain representations of the symmetric group as 
indicated above. The Cj are products of the Clebsch-Gordan coefficients of 5(4). 

C. Formal matrix element calculation 

Given the discussion above, the basic calculation of an interaction matrix element is 
straightforward. We begin with initial and final states defined as 

* = J2^M{a})t,{{T})^,{{r}) r = J2Ck4{W})mT})mr}) (21) 

j k 

We then formally calculate the matrix element to give 
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a</3 k I 



(22) 



j,k ^ 

In writing this, we have integrated formally over spin and isospin coordinates to develop 
the interaction matrix element into simpler spatial matrix elements involving individually 
position-dependent potentials V^'^ given by 



V^'{{r}) = qcJ s,{{a})tj{{T}) 



s'kiWmir}) 



(23) 



In this way we can reduce the spin and isospin dependent nuclear matrix element into a set 
of simpler matrix elements, each one of which involves a position-dependent potential. 
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IV. BASIC VACUUM MATRIX ELEMENT CALCULATION 



Before attempting a calculation of a specific interaction matrix element in the lattice case, 
it is appropriate to first examine the simpler equivalent calculation in the vacuum case. For 
this specific calculation, we will select a four-body problem in which there are two deuterons 
in the initial state, and in which the final state contains a three nucleon body (^He or t) and 
a one nucleon body (n or p). 

A. Vacuum nuclear wavefunctions 

The specific matrix element that we select for our example involves a quintet {S — 
2, Ms — 1) initial state with two deuterons (the deuterons each have spin 1), and a 
singlet {S — 0, Ms — 0) final state which is a hnear combination of + p and ^He + 
n. The reason for this is that a linear combination of states with different total isospin 
{T — and T — 1) are required to resolve a single proton or neutron in the final state 
channel. For simplicity, we have adopted initial and final state wavefunctions with total 
isospin T — 0. As we discussed briefly above, nuclear wavefunctions can be constructed 
from linear combinations of products of basis vectors of representations of the symmetric 
group. For the case of four-body wavefunctions, we have indexed the specific basis vectors 
in terms of their associated Yamanouchi symbol in Appendix B. We have also listed the 
results of a systematic construction of all such four-body wavefunctions with total isospin 
zero in Appendix C. 

The initial state is constructed according to 

*i = *6 = -^[4^5Slot6 - Aswt^] (24) 

where is the notation used in systematic construction of Appendix B. Each of the four- 
particle spin (s), isospin {t), and spatial (ip) functions are individually basis vectors of 
representations of the symmetric group listed in Appendix B. The final state is constructed 
according to 

= ^1 = ^10 Me - Sets] (25) 

where is also the notation used in Appendix B. We have selected these states in particular 
due to the close correspondence between the physical states discussed below, and the group 
theoretical Yamanouchi basis states for these channels. In the case of other channels, a 
linear superposition of Yamanouchi basis states may be required. 
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B. Matrix element in terms of Yamanouchi basis functions 



The matrix element can be evaluated using the approach discussed above. We first express 
the matrix elements in terms of Vhj{1,2), since the Yamanouchi basis functions are either 
symmetric or antisymmetric under 1^2. We write 



6{^f\VHAl,2)\^,) 



(26) 



Using the initial and final states in the form listed above, we obtain four different contribu 
tions to the matrix element 



M 



fi 



(s5^6''/'lo|Vffj(l,2)|siot6'05)-(s5i6^1o|Vffj(l,2)|siot5^6) 



- {s6h^plo\VHJ{l,2)\slotQ^p5) + (sets^/'iol Vhj(1, 2)|siot5V'6) 



(27) 



Because of the symmetry of the interaction and the wavefunctions, only one term survives 
the spin and isospin algebra. The details of this calculation for the tensor term are worked 
out in Appendix [Dl The generalized centripetal term can be worked out in a completely 
analogous way. The net result is that 



Mf^ = (V^lo|V^°''(l,2)|V^5) 



(28) 



where 



V^'^\l,2) = 18v^^^:i^^fl^y-(;.ri,) + ^{L^U + U^)ytL{fin,) (29) 

r^2 IT' 

where and Z/+ are relative 12 (read as one-two) angular momentum operators. 

We have already noted that the Yamanouchi basis is naturally purely symmetric or an- 
tisymmetric under the exchange 1 2 (This is further discussed in Appendix [D]). As a 
consequence, the calculation of the many-body matrix element leads to the simplest al- 
gebraic result when specified in terms of 12 terms. Interactions between all particles are 
represented, since the ipio and ipt^ representations include permutations of the particle num- 
bering. Hence, the use of a position-dependent interaction potential given in terms of 
is appropriate, and implies the use of ri2 as a favored relative coordinate to use in the 
evaluation of the spatial matrix element. 



14 



C. Physical wavefunctions 



The group theoretical machinery that we used in the evaluation above results in a com- 
pact expression given in terms of spatial Yamanouchi basis functions. As such, they have 
certain symmetries with respect to particle interchange. These must be implemented in the 
construction of the spatial wavefunctions. This can be accomplished by using an appropriate 
superposition of permutations of wavefunctions with no particular symmetry. 

For example, suppose we begin with an unsymmetrized "physical" wavefunction for the 
two deuteron initial state 

V^(12; 34) = - ri) 0,(r4 - Vs) F,,, (^4^, ^4^) (30) 

In this wavefunction, nucleons 1 and 2 make up one deuteron given by 0d(r2 — ri), and 
nucleons 3 and 4 make up another deuteron 0d(r4 — fa)- The generalized channel separation 
function F2^2 is the probability amplitude associated with the two center of mass coordinates 
of the deuterons. As discussed in Appendix E, we can use this unsymmetrized physical 
wavefunction as the basis for the construction of the associated Yamanouchi tf]^ function, 
which leads to 



V'5 = 



2iP{l2] 34) + 2^^(34; 21) - -0(23; 14) - -0(14; 23) - -0(13; 24) - -0(24; 13) 



(31) 



This generates the appropriate Yamanouchi basis vector for use in the matrix element ex- 
pression above. 

In the case of the final state channel, we begin with the physical wavefunction 



V^(123; 4) = 03//e(r3 - r2, r2 - ri) Fg^i 



ri r2 + ra 



(32) 



where 03i7e is the ^He wavefunction specified in terms of two relative internal coordinates, 
and where F^^i is the generahzed channel function for the center of mass for the ^He and 
neutron. Prom this physical wavefunction, we can construct the associated Yamanouchi ■0io 
function according to 
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^^(123; 4) + V'(124; 3) + ^'(134; 2) + ^'(234; 1) 



(33) 



This function can be used in the group theoretical matrix element formula above. 
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V. LATTICE RESULT 



Let us now consider exactly the same physical process taking place inside the lattice. 
We are interested in the calculation of the interaction matrix element with wavefunctions 
that contain both nuclear and lattice coordinates. Perhaps the simplest way to approach 
the problem is to make use of the vacuum results from the last section, but expand the 
definition of the initial and final states to include the lattice. 



A. Wavefunctions for lattice and nuclei 

We now augment the definition of the physical wavefunctions to include the lattice coor- 
dinates. For example, in the new physical initial state wavefunction, we include the center 
of mass coordinates of the other nuclei to obtain 



V'(12;34,{R}) = 0.(r2-rO0,(r4-r3)F^,(^i±^,^i±^,{R}) 



(34) 



The other nuclei in the lattice are spectators in the sense that they do not participate in 
nucleon rearrangement associated with the reaction. Consequently, the appearance of the 
associated degrees of freedom into the lattice channel separation factor is an appropriate 
generalization of the generalized channel separation factor from the vacuum case. We still 
require the nucleons that participate in the reaction to be described by a wavefunction that 
is totally antisymmetric under the exchange of nucleons. The vacuum formulation that we 
discussed above accomphshed this in the absence of spectator nuclei. There seems to be no 
reason that we cannot adopt precisely the same kind of formulation here (this is discussed 
in Appendix F). Hence, we form the appropriate Yamanouchi spatial basis function for the 
initial state made up of the same superposition of physical states that are augmented with 
spectator coordinates. This leads to 



^5^ 



2V^(12; 34, {R})+2^(34; 21, {R})-^(23; 14, {R}) 



^{U; 23, {R}) - ^^(13; 24, {R}) - ^(24; 13, {R}) 



(35) 



Final state wavefunctions can be developed similarly. The physical wavefunction for the 
final state is 

V^(123; 4, {R}) = 03He(r3 - r„ r, - r,) F^^, (^^1±^±^, r,, {R}^ (36) 
The appropriate final state Yamanouchi basis is formed through 



16 



^10 



V^(123; 4, {R}) + V^(124; 3, {R}) + ^^(134; 2, {R}) + ^(234; 1, {R}) 



(37) 



B. Matrix element 

Wc can take advantage of the same group theoretical construction to evaluate the matrix 
element when spectator nuclei are present, since we only required total antisymmetrization 
on the four nucleons involved in the local strong force interaction. Consequently, we may 
write 



(38) 



We see in this result a very strong connection between the vacuum result and the lattice 
result. In general we can systematically extend vacuum calculations to the lattice case 
through the addition of the lattice nuclear center of mass coordinates as spectator degrees 
of freedom. 



C. Matrix element in terms of physical wavefunctions 

We have made use of representations of the symmetric group in order to construct nuclear 
wavefunctions and to evaluate the interaction matrix element. In the case of the spatial 
wavefunctions, these basis functions are made up of several different permutations of the 
physical wavef unction. If we make use of these expansions, we can recast the interaction 
matrix element in terms of matrix elements involving physical wavefunctions. This will in 
general produce a large number of terms. 

In the specific case of the tensor matrix element that we have been using as an example, 
we obtain twenty four terms in total. We may write 

V^'^{r2-ri)(t)d{v2-ri) <pd{rA 

(^4^, ^4^, {R}) d^v, d'v, d^ra d^v^d^^'^B} + • • • (39) 
where the • • • represent the other twenty three terms in the expansion (see Appendix G). 



V12I 



hHe{rz - r2,r2 - ri) F; 



L / ri + r2 + ra 



3,1 



,r4,{R} 
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VI. MATRIX ELEMENT IN TERMS OF PHONON COORDINATES 



A primary goal of this paper is to develop an expression for one of the integrals that 
appear in the matrix element explicitly in terms of phonon coordinates, similar to what we 
discussed in Section II. We begin by defining the integral / according to 



-I 



^2,He(yz - r2,r2 - Ti) F^i 



ri + r2 + rs 



,r4,{R} 



yR'^{^2 - ri)^/)d(r2 - ri) (?!)d(r4 - r^) 

(40) 



A. Integral in terms of center of mass and relative coordinates 

This integral includes integrations over nucleon coordinates as well as the center of mass 
coordinates of the spectator nuclei. We would like to recast the integral in terms of center of 
mass coordinates and relative coordinates, to set things up for a transformation to phonon 
mode coordinates. We may write 

d'^i2 dHLd^^He dH dHil d^is d'Rii, d^R^i^ d'^'-^R} (41) 
where we have defined the relative coordinates 



I21 



^43 



and the center of mass coordinates 



R 



r\ + ri 



■dl 



R 



S32 



4 + 4 



■d2 



R^ 



S21 



f f f 

r{ + + 



(42) 



(43) 



B. Deuteron-deuteron separation 

When the two deuterons are close enough for strong force interactions to take place, then 
the Coulomb repulsion is sufficiently strong that this interaction dominates over the forces 
from other atoms. Although we can use a phonon description in this case, it would not 
be the most natural, and the Coulomb forces would require the presence of very high-order 
phonon operators. Hence, it may be more convenient to adopt the point of view that the 
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two deuterons appear to the lattice as an equivalent He nucleus when they are a few fermis 
apart. The separation between them may best be handled then as part of the microscopic 
nuclear problem. In this case, we may write 

FI2 - RW + {R}) 5^^'^ (^-1^2 - l^L + Rk - (-^^21 + Rk - \^m)) 

S^'^ [-1^2 + l^L + Rk - (^^21 + Rk - Ir^)) 

S^''^ (^^32 + + Rk - (-^^43 + RW + ^r^.)) 

S^'^ (r{ - (^^^3 + RWe + ^r^dd)) d'^i2 d'^id'^ine d^ d'^l, d^Rl^e dh^,,d'''-\K} 

(44) 

In writing this, we have made use of the initial state relative and center of mass coordinates 

R-We = 2 ("^'^i "^'^2) ^dd = R-d2 - R-dl (45) 

Note that in casting the matrix element in this way we are not assuming that the mi- 
croscopic ^He nuclear wavefunction is the same as the microscopic nuclear wavefunction for 
two deuterons. Instead, we recognize in these expressions that there is little difference from 
the point of view of the rest of the lattice between two deuterons localized on the fermi 
scale (which the rest of the lattice sees as a charge 2 and mass 4 object), and a ^He nucleus 
(which the rest of the lattice also sees as a charge 2 and mass 4 object). Due to the Coulomb 
repulsion between the two deuterons when they are close, the associated relative wavefunc- 
tion is far into a tunneling regime. We have the choice to describe this tunneling in terms 
of deuteron coordinates as phonon operators, which is complicated, or else to describe this 
tunneling via a relative nuclear coordinate based on a ^He center of mass coordinate, which 
is much simpler. The formal problem is not changed with this replacement, although the 
approximations used in the evaluation of the resulting expressions could be different. 

C. Transformation to phonon coordinates 

We are now in a position to make the transformation to phonon coordinates. In the initial 
state, a nuclear center of mass coordinates can be expressed in terms of phonon coordinates 
in general according to 

= + E"}N?;. (46) 

m 



/ <p3Hei^i2, ^21) ^3^:1 (Rk' r{, {R}) 



19 



Consequently, we may transform from initial state center of mass coordinates to initial state 
phonon coordinates 



(47) 



where is the center of mass of the initial state lattice. A similar transformation is possible 
for the final state center of mass coordinates, which we indicate by 



where R{ is the center of mass of the final state lattice. 



where 



- ^rj,„ + ^rj,„ {R^}^ 6^'^ (r{ - Rl) 

S^'^ {^1^2 + + ^Le - il^ll + RW - 

S^'^ (^^2 + + Rk - (-^^^3 + Rii/e + ^rj,,)) 



d'Ri d'B!^d'e32 d'e2i d'ri d"^, d'e.s d^r^^d'^^-^fd 



'3Af-3qi 



(48) 



(49) 



R^ 

H 

R" i 
me 



R{ (q^R{) 
RW Rl) 

Ri = RUq\R';,) 



(50) 

(51) 

(52) 
(53) 



D. Integral in terms of generalized lattice wavefunctions 

The generalized lattice channel separation factors are the same as lattice wavefunctions 
augmented with additional degrees of freedom. Let us adopt a notation of the form 
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2,2 



R 



^/(q^r{) 



(54) 



(55) 



The generalized lattice wavefunctions include an explicit dependence on phonon coordinates, 
as well as an explicit dependence on a microscopic nuclear degree of freedom. In the initial 
state, this extra degree of freedom is the deuteron-deuteron separation; in the final state it 
is the position of the neutron. Implicit in this description is an assumption that the lattice 
center of mass is preserved, which is an additional assumption here that is consistent with 
our focus on phonon exchange (we do not mean to imply here that no coupling is possible 
with the lattice center of mass coordinate). With this notation, we may write 



*^(q/,r{)j%;,(q/,q^)vI/f(q\rj,,)5(37V-3)^„/ 



A • q^ - h)d^ri d^r^rfdq-^dq* 



(56) 



where the interaction potential is 



*/i(q^,q* 



e 2 dd' 



^(3) 



ri - (^^43 + ^He + ^r^dd) ) ^'^^2 d'$L dHi, 



We may integrate over the final state phonon coordinates to obtain 



7iq^r{) 



Vfiiqf', q')^i:{q\ rj,Jrf'r{ dh^iJq 



(57) 



(58) 



in which the final state phonon coordinates are expressed in terms of initial state coordinates 



as 



Aq* 



(59) 



Equations ( 1571) and (l58l) are the primary results that we find in the particular integral 
involved in the specific deuteron-deuteron matrix element considered as our example. Similar 
results can be readily obtained for other integrals appearing in this matrix element, and also 
for other matrix elements as well. 
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VII. CONCLUSIONS 



The analysis presented in this paper is motivated by claims of excess power production 
in experiments involving metal deuterides that have been reported over the years since the 
original claim of Fleischmann and Pons in 1989. The existence of such an effect implies new 
physical mechanisms, in particular mechanisms which lead to a coupling between condensed 
matter degrees of freedom and nuclear degrees of freedom. There has been confusion and 
conjecture about how energy or angular momentum exchange with the lattice works, and 
no relevant formulation or calculations have been reported in the mainstream literature. 
In recent years, our focus has been on a particular scheme involving excitation transfer, 
idealized models for which are discussed in detail in a recent manuscript.- The formulation 
discussed in this paper leads to formulas that can be used for a quantitative evaluation of 
this model, and other models that have been proposed which also involve coupling with 
phonons. 

The nuclear part of the calculation requires spin and isospin degrees of freedom in addition 
to the spatial degrees of freedom. We developed nuclear wavefunctions based on Yamanouchi 
basis functions for each degree of freedom individually, which allows us to integrate out 
the spin and isospin degrees of freedom separately in a convenient way. This leads to a 
straightforward expansion of the interaction matrix element in terms of a sum of spatial 
integrals which contain different permutations of particle numbering. 

The lattice part of the problem comes into the calculation initially as spectator degrees 
of freedom in the generalized channel separation factor. In essence, the microscopic part of 
the calculation at this level is completely equivalent to the vacuum calculation, since the 
spectator nuclei are not involved locally in the strong force interaction. The lattice does come 
into the problem in a more fundamental way when we replace the nuclear center of mass 
degrees of freedom with phonon coordinates. Since the reacting nuclei in this formulation 
are taken to be part of the lattice in both initial and final states, the lattice is changed as 
a result of the reaction. We chose a description in which both the initial and final state 
product nuclei are assumed to be included in the associated lattice definitions. This leads 
to integrals which contain both initial state and final state lattice wavefunctions, where the 
lattice structure is different in the two cases. 
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APPENDIX A: THE NUCLEAR POTENTIAL 



We make use of the Hamada- Johnston potential as our model nuclear potential in this 
paper. It is one of the simplest isospin-preserving nuclear potentials to give reasonable 
results. The Hamada- Johnston potential between nucleon 1 and nucleon 2 can be written 

as 



Vhj{1,2) = Vc{l,2) + Vt{1,2)S,2 + Vls{1,2){L ■ S) + Vll{1,2)L,2 (Al) 



where L and S are the relative angular momentum and and the total spin of the two nuclei. 
The individual potentials are defined as 



Vc{l, 2) = (ti • T2) {<Ti ■ 0-2) yc{l^ri2) 
Vt{1, 2) = (ti • T2) yT(l^ri2) 



Vls{1,2) = |/i5(m2)(L-S) 

L12 = (o-i • (T2)V - -{(Ti ■ L)(<T2 • L) - -(<T2 • L)(<Ti • L) 



where 171.,^ is the pion mass, a stands for odd-singlet (os) , even-singlet (es), odd-triplet (ot) 
or even-triplet (et). The functions Y{x) and Z{x) have the definitions 



The spatial functions are given by 




Y{x) 



e 



X 



—X 
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APPENDIX B: SYMMETRIC GROUP 



The symmetric group S{n), the group of permutations of n objects, plays an important 
role in our calculation. It does so, in three important ways 

• Clebsch-Gordan coefficients are used to construct the antisymmetric wavefunctions. 

• Schur-Weyl duality lets us build wavefunctions with a well-defined spin (and isospin) 
and which also transform as Yamanouchi basis vectors under S{4). 

• Induction coefficients of S'(4) are utilized to construct Yamanouchi basis vectors for 
the spatial part. 

In quantum mechanics, we are familiar with the transformation properties of the spherical 
harmonics. In the case of Yi^s the / labels the irreducible representation of 5*0(3) and m 
labels the basis vector. Similarly for the symmetric group, we get two labels: the Young 
diagrams and the Yamanouchi symbols22i2i^ where the Young diagrams label the irreducible 
representation of 5'(4) and the Yamanouchi symbols label the basis vectors. The index 
scheme we use for the Yamanouchi symbols is defined in Table [B 



TABLE I: Definition of indices for the Yamanouchi symbols 



Yamanouchi Symbol 


index 


Yamanouchi symbol 


index 


4321 


1 


2121 


6 


3211 


2 


2111 


7 


3121 


3 


1211 


8 


1321 


4 


1121 


9 


2211 


5 


1111 


10 



1. Clebsch-Gordan coefficients 

The general formula for constructing antisymmetric wavefunctions is given by 

^imnif"]) = E ^Iriv'my'' s{[nY')T{[nY") (bi) 

Y,Y ' ' ' Y',Y" 

where the C's are the Clebsch-Gordan coefficients of 5'(4), [/]'s are the Young diagrams, F's 
are the Yamanouchi symbols, / and Y are the conjugate Young diagrams and Yamanouchi 
symbols and [1^] represents the completely antisymmetric one dimensional representation 
with one Yamanouchi basis vector (index 1). The Clebsch-Gordan coefficients of 5'(4) are 
available in the Iiterature^>2ii22ii2^. 
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As a specific example, we construct the basis function = "^q. Using our indexing 
sclieme of Table [T] above and Table 4.13 of Chen et aP^, we see that 



Cio,5 = 1 C'lce = 1 (B2) 
Hence the spin-isospin part of Equation IBll yields 

{St)5 = Slots (-5^)6 = SiotQ (B3) 

where the subscripts in (5^)5 or {st)^ again refer to the Yamanouchi basis for the particular 
spin-isospin combination. In addition, we have 

C,---L (B4) 

Therefore, the totally antisymmetric wavefunction can be written as 

" "^^5(st)6 - ;^V^6(st)5 = [^5^10^6 - ^6^10^5] (B5) 



2. Schur-Weyl duality 

The Schur-Weyl duality relates certain representations of SU{m) or GL{m) to the rep- 
resentations of S{n). For our purposes, since we are only interested in the spin-half (or 
isospin-half) case, the Schur-Weyl duality is particularly simple. It states that there are no 
spin (or isospin) wavefunctions corresponding to the index 1,2,3 or 4. It also allows us to 
construct the rest of the wavefunctions very simply. In the usual spectroscopic notation, 
with the nucleon number as a subscript and a referring to spin angular momentum. 



S5 = 




= l),0'3('S'l23 = 


5) 


, 0-4 


s = 


0) 


s& = 


0-10-2(5*12 


= 0), 0-3(^123 = 


^) 


, 0-4 


s = 


0) 


S7 = 


0"lO"2(5'i2 


= 1), 0-3(5'l23 = 


5) 


, 0-4 


s = 


1) 


58 = 


aia2(S'i2 


= l),0-3(5'l23 = 


^) 


, 0-4 


s = 


1) 


S9 = 




= 0), 0-3(^123 = 


5) 


, 0-4 


s = 


1) 


SlO = 


0-10-2(5*12 


= l),0"3(5l23 = 


|) 


, 0-4 


s = 


2) 



Here on the RHS, we are suppressing the value of Ms because different Ms values transform 
in the same way under the symmetric group. Exactly the same results apply to the isospin 
wavefunctions, since they are also angular momentum 1/2 representations of SU{2). 
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3. Induction Coefficients 



In the previous subsection we explicitly constructed the spin/isospin Yamanouchi wave- 
functions using angular momentum addition. We still have to construct the spatial Ya- 
manouchi wavef unctions. One way to do it to use induction coefficients of the symmetric 
group, which are tabulated and discussed in Chen et al.— . The relevant formula for our 
purposes is 

riUlYf) = Yl Sll]'^iMy'^-,.r\[9],Y^^ (B7) 

Here n = nl +n2, [f] and are the Young diagram and the Yamanouchi symbol for an n- 
particle spatial wavefunction, [g] and are the Young diagram and Yamanouchi symbol for 
the nl-particle wavefunction and [h] and Y^ are the Young diagram and Yamanouchi symbol 
for the r;,2-particle wavefunction and u = {uji, 002) represents a permutation of {1, 2, ■ ■ ■ , n}. 

As a concrete example let us consider n = A case with nl = n2 = 2. Suppose we want 
to construct a ipQ, starting with two two-body spatial wavef unctions, and 0', which are 
completely symmetric. Hence the summation over the Yamanouchi symbols is trivial and 
we only need to sum over various tu's corresponding to different orderings of {1, 2, 3, 4}. We 
can use Table 4.17 of Chen et aF^ to see that (suppressing the trivial Yamanouchi symbols) 

1 111 

06 06 q6 q6 ('RC'\ 

•^{2,3,1,4} ~ 2 {1.4,2,3} — 2 '-^{1,3,2,4} " 2 '^{2,4,1,3} " 2 ^ ^ 

Hence 

7/>6 = ^ [-0(2, 3)0'(1, 4) - 0(1, 4)0'(2, 3) + 0(1, 3)0'(2, 4) + 0(2, 4)0'(1, 3)] (B9) 

We can use this approach to develop a complete set of results needed for the case of two 
deuterons as used in this paper. For the initial state, we may write 

V's = ^ [2V'(12; 34) + 2^^(34; 21) - ^(23; 14) - V(14; 23) - V(13; 24) - ^^(24; 13)] 
V 12 

= ^ [-0(23; 14) - 0(14; 23) + V(13; 24) + V(24; 13)] (BIO) 

where "0(12; 34) is any spatial wavefunction which is symmetric under the 1^2 and 3^4 
exchange. For the final state, we are using a three-nucleon body and a one nucleon body. 
The Yamanouchi wavef unctions are 
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V'lo = ^[^(123;4)+^(124;3)+V(134;2)+V(234;1)] 
V'r = 4^[3^(123;4)-V(124;3)-V(134;2)-V'(234;1)] 



V's = -^[2V'(124;3)-V(134;2)-^(234;1)] 

V^g = ^ [^-(134; 2) - V(234; 1)] (Bll) 

where the -0(123; 4) is any wavefunction symmetric under any permutation of {1, 2, 3}. 

Hence we see that by using induction coefficients, we can take spatial wavefunctions 
which belongs to certain representations of subgroups of 5'(4) and use them to construct 
Yamanouchi basis vectors e.g. for the case of two deuterons, we started with -0(12; 34) which 
is symmetric under the exchange of 1 <-> 2 or 3 4, and constructed spatial wavefunctions 
with Yamanouchi basis index 5 or 6. 
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APPENDIX C: YAMANOUCHI BASIS STATES 



Below we explicitly write down the T = basis states. Please note that the subscripts on 
the LHS of the equations are merely ways of enumerating the possible wavefunctions, where 
as the subscripts on the RHS are the indices for the Yamanouchi symbols (see Table [T]). 
S = 

^1 = ^w-^ Me - Sets] 



^2 = V^i-^ [s5h + Sate] 

^3 = -^^5 [^5^6 + Sets] 



[sbh - sete 



S = 1 



V6 



F^8 



V2s 

7^6 + Sste + Sgts 



2^/3 
1 

1 

2^3 



^2 



2^3 

— Ai/2s7t5 + Sgts — Sgt 



+ 



S = 2 



V2 
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APPENDIX D: THE TENSOR MATRIX ELEMENT 



We have calculated a complete set of such Yamanouchi matrix elements for the T = 
case. There are too many cases to be presented here. We therefore focus on a specific 
example by calculating the tensor force matrix element of \Efi(Ms = 0) and "^qIM^ = 1). 
The main issue with all these calculations is that the parameters in most of the nuclear 
potentials are fit separately to spin singlet-triplet and even-odd spatial parity (for example, 
in the Hamada- Johnston potential the y functions are parameterized by a and b parameters, 
which depend on the spin and parity). The most convenient way of dealing with these issues 
is to use projection operators to separate the various parity and spin channels. However we 
can, without the use of projection operators, directly exploit the properties of the interaction 
and the wavefunctions to evaluate these matrix elements. 



1. Exchange properties of the interaction and wavefunctions 

First of all, since we are using completely antisymmetric wavefunctions, the tensor part 
of the matrix element calculation simplifies to 



y^^VT{a,P)Sap 

a<f3 



VI/,: ) = {^mi Yl VT{^,J)S.,\^e{l)) 

l<j<j<4 

= 6(^i(0)|yT(l, 2)^121^6(1)) 



(Dl) 



Here the values and 1 in the parentheses represent the Ms values of the initial and final 
wavefunctions. So now, as discussed in the main text, the tensor interaction can be evaluated 
by 



VTia,P)Sal3 

a<l3 



{s5t6lplo\VT{l, 2)5i2|siot6V'5)-(s5i6V'lo|W'(l, 2)5'i2|siot5?/'6) 



(s6t5V'io|1^T(l,2)5i2|siot6'05) + (sgtsV'io I Vt(1, 2)5i2 IsiotsV'e) (D2) 



From the explicit form of Vt(1, 2)5*12, we can see that it is symmetric under the exchange 
1 <-i> 2, individually for spin, isospin and space parts. If we look at the explicit form of 
the spin wavefunctions in Appendix [B] we can see that since particles 1 and 2 couple to 
form a triplet, S5, sy, Ss and sio are even under 1^2. Similarly since particles 1 and 2 
couple to form a singlet, sq and sg are odd under 1^2. The same applies to the isospin 
wavefunctions. Being even or odd under the exchange 1 2 in this case has nothing to do 
with angular momentum algebra, but is a property of the Yamanouchi symbol. Hence the 
spatial wavefunctions ip^, ipj, ips and ipio are even under 1^2 and ipQ and ipg are odd under 
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1^2. As a concrete example, this can easily be seen by explicitly looking at the ipQ from 
Appendix [B] 



^6=2 [-</>(2, 3)0'(1, 4) - 0(1, 4)</)'(2, 3) + 0(1, 3)0'(2, 4) + 0(2, 4)0'(1, 3)] (D3) 
This ipQ, under the exchange 1^2, transforms to 

i^e^l [-0(1, 3)0'(2, 4) - 0(2, 4)0'(1, 3) + 0(2, 3)0'(1, 4) + 0(1, 4)0'(2, 3)] = -^g (D4) 

Because of the symmetry of the interaction under 1^2, individually for space, spin and 
isospin parts, we only get non-zero couplings when the space, spin and isospin have the same 
parity (under 1^2) for both the initial and the final wavefunction. Hence we can see that 

(s5i6^1o|^T(l,2)S'i2|siot5^6) = 
(s6)f:5^1o|V'T(l,2)S'i2|Siot6^5) = 
{s6t5^1o\VTil,2)Si2\Siot5^6) = 

The only non-zero matrix element is (ssteV'iol Vt(1, 2)S'i2|siot6'05)- 



2. Matrix element calculation 

We have reduced our matrix element calculation to just calculating 
(s5t6'0io|^(l? 2)S'i2|siot6'^5)- We can evaluate the spin and isospin pieces separately 



{t(i\Ti.T2\tG) 



(D5) 



\S12\s10iMs = 1)) = -2v^zi2 



xi2 + iyi2 



(D6) 
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Since the Sio and S5 are triplets and ipiQ and -05 are even, the ?/t(/^'^i2) should be the one 
parameterized by even parity and triplet spin. Hence yrifi'r 12) = 1/t(/^^i2) where et stands 
for "even, triplet". 

We can assemble these results to obtain 



^i(O) 



l<i<j<4 



^-6(1) ) = 18^3 / ^ 



^10 



{xi2 + iyi2) Z12 VTifJ-ru) 



12 



■05 d^ri ■ ■ ■ (fr^ 
(D7) 

where ipio = V'io(''i; ' ' ' > ''4) V's = V'slri,-'' y'^'i) are functions of spatial coordinates 
alone. 
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APPENDIX E: VACUUM CASE ANTISYMMETRIZATION 

In Appendix [D] we saw that it was straightforward to calculate the matrix element of 
the tensor force in terms of the Yamanouchi basis. This is in general true for all nuclear 
interactions. Perhaps the simplest way of carrying out realistic physical calculations is to 
express initial and final state wavefunctions in terms of the Yamanouchi basis. Then we can 
use the Yamanouchi basis matrix element results to calculate the nuclear matrix elements 
with physical initial and final states. 

For an isospin preserving, T = reaction, in which two deuterons interact to form a 
three-nucleon body and a one nucleon body the relevant states are: 

• ^2,2: One initial 2 + 2 state of two deuterons in a quintet spin 2 state. We should 
in general consider the deuterons to be in a singlet or triplet as well. However, we 
are focusing on a single example here. In this enumeration of states, we suppress the 
various Ms values. 

• \I^3,i: Two final 3 + 1 states (since these can be singlets or triplets). These are a linear 
superpositions of the + p and ^He + n states. 

1. Physical wavefunctions 

The wavefunctions discussed above are formal Yamanouchi objects which are very general. 
We need to focus on specific wavefunctions for the different mass 4 (T = 0) channels in order 
to proceed. For simplicity, we adopt a wavefunction for two deuterons that are in a quintet 
(spin 2) state. In this case, we can use 



Here, A is an antisymmetrizer, and ip{12\ 34) is the spatial part of the deuteron wavefunc- 
tions. This wavefunction can be taken to be of the form 



There are two kinds of 3 + 1 wavefunctions, including singlet and triplet states. The 
singlet 5* = states can be written as 



^2,2 



^{^(12;34)siot6} 



(El) 




(E2) 




(E3) 



The 3+1 triplet 5* = 1 states can be written as 




(E4) 
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The spatial part of ^3,1 and ^3 ^ is of the form 



V'(123; 4) = 03i/e(r3 - r2, - ri) ^3,1 I ''^^'L'^''' , r4 ) (E5) 



2. Physical wavefunctions in terms of Yamanouchi basis states 

All the physical wavefunctions for the initial and final states can be expressed as hnear 
combinations of the Yamanouchi functions \l/4 and ^g- This is so because of group 
theoretical considerations and our assumptions that the spatial part of the deuteron, triton 
and helium wavefunctions are symmetric under the exchange of any two particles. The 
results in particular are 

*2,2 = *6 (E6) 

with 

V's = [2^(12; 34) + 2^^(34; 21) - ^(23; 14) - ^{U; 23) - ^(13; 24) - ^(24; 13)] (E7) 



^6 = 1 [-V'(23; 14) - ^^(14; 23) + V(13; 24) + V(24; 13)] (E8) 



In the case of the singlet 3+1 channel, the antisymmetrizer acts to produce a single Ya- 
manouchi basis state 

^3,1 = ^1 (E9) 
where the associated fully symmetric spatial part ipio can be expressed as 

^10 = ^[V^(123;4) + V'(124;3) + V'(134;2) + ^(234;l)] (ElO) 

In the case of the triplet 3+1 channel, the antisymmetrizer acts to produce a single Ya- 
manouchi basis state 

^3,1 = ^4 (Ell) 
where the associated mixed symmetry spatial parts ip^, ipg, and ipg are 

^7 = ^ [3^(123; 4) - V(124; 3) - ^(134; 2) - ^(234; 1)] (E12) 
V 1^ 



= ^ [2V'(124; 3) - ^^(134; 2) - ^^(234; 1)] (E13) 
= ^ mm 2) - V(234; 1)] (E14) 
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APPENDIX F: LATTICE CASE ANTISYMMETRIZATION 



In this appendix we consider a partial antisymmetrization in the case of lattice wavefunc- 
tions. We retain the basic vacuum description for the four nucleons of the vacuum case, and 
extend the description to include the other nuclei as spectators. 



1. Physical Wavefunctions 



We can implement this program in the case of the initial states to write 

*2,2 = ^{^^(12;34,{R})siot6} (Fl) 

Here, A is the antisymmetrizer acting on particles 1, 2, 3 and 4 only, and ip{12; 34, {R}) is 
the spatial part of the deuteron and the lattice wavefunctions. This wavefunction can be 
taken to be of the form 



^(12; 34) = 0,(r2-rO0.(r4-r3)Fi;2(^^li^,^l±^,{R}) 



(F2) 



There are two kinds of 3 + 1 wavefunctions, including singlet and triplet states. The 
singlet S — states can be written as 

= ^ {V'(123; 4, {R})^{s,te - s^t,)"^ (F3) 
The 3+1 triplet S — 1 states can be written as 

%^ = A 1^(123; 4, {RD^Me - ^9^5)} (F4) 
The spatial part of ^3 ^ and ^'^i is of the form 

^^(123; 4, {R}) = cl>sHe{rs - v^, - v,) F^, f^l±^±il, r4, {R}) (F5) 



2. Physical wavefunctions in terms of Yamanouchi basis states 

Now since our antisymmetrizer, acts only on the nuclear coordinates, and not on the 
{R}, we see that the results of the vacuum case hold except that the spatial parts of the 
wavefunctions get augmented. For example, we may write 

*2,2 - *6 (F6) 

with 
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V's = [2VX12; 34, {R}) + 2V^(34; 21, {R}) - V(23; 14, {R}) - V(14; 23, {R})- (F7) 



V^(13;24,{R}) -V(24;13,{R})] (F8) 
V'6 = ^ [-^-(23; 14, {R}) - V(14; 23, {R}) + ^(13; 24, {R}) + V'(24; 13, {R})] (F9) 

In the case of the singlet 3+1 channel, we obtain a similar connection to the Yamanouchi 
basis as in the vacuum case 

= *i (FIO) 
where the spatial part is now augmented to include the spectator nuclei 

V'lo = I mi23; 4, {R}) + ^^(124; 3, {R}) + ^(134; 2, {R}) + ^^(234; 1, {R})] (Fll) 
The triplet 3+1 channel becomes 

= *4 (F12) 

where 



^7 = [3V'(123; 4, {R}) - ^(124; 3, {R}) - ^(134; 2, {R}) - ^(234; 1, {R})] (F13) 

viz 



V-s = ^ [2V(124; 3, {R}) - ^(134; 2, {R}) - ^^(234; 1, {R})] (F14) 
V-g = ^ [V'(134; 2, {R}) - V(234; 1, {R})] (F15) 
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APPENDIX G: INTEGRALS RESULTING FROM PERMUTATIONS 



Due to the antisymmetrization, we generate a total of twenty four terms in the matrix 
element example in the main body of the paper. These are 



12 (*3,i(M. = 0)|F(1,2)|*2,2(M, = 1)) = 

J V(123;4rt;(l, 2)^(12; 34)+ J V(123; 4)* ^1, 2) - 1^1, 4) - ^t;(4, 1)]V(12; 43) 
+ J V(123;4)*[t;(l,4)-i;(l,2)]V(14;23)+ J ^(123; 4)*i;(4, 1)^^(41; 23) + J V(123; 4)*t;(l, 2)V'(34; 12) 
+ J V(123;4)*Kl,2)-l^;(l,4)-i^;(4,l)]V^(43;12)+ J V(123; 4)* [-^;(1, 2) + t;(l, 4)]V(23; 14) 
+ J V(123; 4)*i;(4, 1)^^(23; 41) + J ^(123; 4)* [-i;(l, 2)]V'(24; 13) 

13)+ / V'(123;4)*[-i;(l,2)]V(13;24) 



+ J V(123;4)*[-^^;(l,4)-^t;(4,l)]V(42 
+ J V'(123;4)*[-li;(l,4)-i^;(4,l)]V'(13 



42) (Gl) 



where V{1,2) is some general nuclear potential which contains spin and isospin dependent 
terms, and where v{i,j) is the spatial interaction between particles i and j which is obtained 
by carrying out the spin and isospin algebra. In the special case that 1^(1,2) = ¥^7(1, 2), 
thenv{i,j)^V^''\i,j). 

Some of these terms are amenable to a simple interpretation in terms of stripping re- 
actions. Others can be considered as higher-order exchange terms. However, in nuclear 
reactions, it is not necessary that higher order exchange terms are necessarily less important 
than simple direct or exchange interactions. Hence, in any realistic computation we will 
need to sum over all of these terms. 
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